Binary Disruption by Massive Black Holes: Hypervelocity Stars, S Stars, and 

Tidal Disruption Events 

Benjamin C. Bromley 

Department of Physics & Astronomy, University of Utah, 115 S 1400 E, Rm 201, Salt Lake City, 
O ! UT 84112 

Hi bromley@physics.utah.edu 

Scott J. Kenyon 

On. 

CN . skenyon@cfa.harvard.edu 

^ • Margaret J. Geller 

o' 



43 

Ok 

6 

CO 



mgeller@cf a . harvard . edu 

Warren R. Brown 

Smithsonian Astrophysical Observatory, 60 Garden St., Cambridge, MA 02138 

wbrown@cf a . harvard . edu 



> ■ ABSTRACT 

OO 

\^q \ We examine whether disrupted binary stars can fuel black hole growth. In this 

^ ' mechanism, tidal disruption produces a single hypervelocity star (HVS) ejected at high 

velocity and a former companion star bound to the black hole. After a cluster of bound 
stars forms, orbital diffusion allows the black hole to accrete stars by tidal disruption 
at a rate comparable to the capture rate. In the Milky Way, HVSs and the S star 
cluster imply similar rates of 10 -5 — 1CP 3 yr _1 for binary disruption. These rates are 
consistent with estimates for the tidal disruption rate in nearby galaxies and imply 
significant black hole growth from disrupted binaries on 10 Gyr time scales. 

Subject headings: Galaxy: kinematics and dynamics — Galaxy: halo — Galaxy: stellar 
content — stars: early-type 



1. Introduction 

Massive black holes with masses M. < 10 9 M inhabit the centers of many galaxies. It 
is uncertain how these obj ects form and grow. Proposed ideas include the direct collapse of a 



primordial gaseous cloud (JEisenhauer et al.1 120051 ; iBegelman et al.l 120061 ) , accretion of gas from 



a surrounding disk (IDebuhr et al.l I201Q ) , capture of stellar- mass objects from a nucle ar cluster 
(JMerritt &: Poonl |2004j) , and hierarchical merging of black holes (jVolonteri et al.l 120031 ) . For all 
mechanisms, growing the most massive black holes is challenging. Gaseous accretion flows are 
limited by the need to shed angular momentum for gas to fall efficiently onto the black hole. Pairs 
of black holes tend to eject material from their vicinity, slowing orbit contraction and sta lling the 



growth of either black hole (JMiralda-Escude fc Gould 



2000 



Milosavljevic fc Merrittl 120031 ) . 



Recent observatio ns of candidate stellar tidal disruption events (TDEs; e.g.. lGezari et al.ll2009l ; 
van Velzen et al.l l201ll ) have renewed interest in stellar capture hypotheses. When a single star 
wanders inside the black hole's tidal radius, t he black ho le shreds the star. Some material ac- 
cretes onto the black hole and powers a flare (IReeslll988r). Flare evolution depends on proper 



ties of the black hole and the shredded star (e.g., lLoeb k, Ulmerlll997l : IStrubbe fc Quataert 



2009 



Lodato h Rossill201ll . and references therein) . Obs ervational and theoretical estimates suggest cap 



ture rates of 10 5 — 10 3 yr 1 per nearby galaxy ( Wang &: Merrittl I2004J : Ivan Velzen et al 



sufficient to grow a fai rly massive black hole on cosmological time scales (e.g.. lMerritt &: Poon 



2011 



2004 



Brockamp et al. 



2011 



Here, we assess the contribution of bin ary stars to TDEs and the growth of bl ack holes. Binary 

1991). When a binary 



stars are a significant fraction of all stars (jAbtl Il983l : iDuquennoy &: Mayor 



system wanders too close to a black hole, the three-body i nteraction e jects one binary partner at 
high speeds and leaves the other bound to the black hole (JHillslll988l ). In addition to predicting 
the ejected binary components as hypervelocity stars (HVSs), the Hills mechanism yields bound 
stars which serve as a mass reservoir for TDEs and black hole growth. Measurements of both the 
bound stars and the HVSs provide direct constraints on the capture frequency along with robust 
predictions for the frequency of TDEs and the rate of black hole growth in the Galactic Center. 

2005, 



We show that this binary model is consistent with observations of HV Ss (Brown et al. 



2012 ). and the S star cluster in the Galactic Center (JGould &: Quillen 



2003 



O'Leary fc Loeb 



2008) 



(§2). We use a simple two-body relaxation model to show how captured stars can diffuse onto orbits 
that intersect the black hole, leading to a rate for TDEs similar to the observations (§3). Finally, 
we demonstrate the potential for the binary capture mechanism to grow supermassive black holes 
(§4). 



2. Binary Disruption Rates from HVSs and S Stars 

To motivate the binary star model, we consider the collision cross-section of the black hole, 
a m . Black holes shred a single star with mass m and radius r at the tidal radius, 
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When M, < 10 8 M , this interaction occurs outsid e the Schwarzschild radius; the black hol e 
accretes a fraction / acc ~ 0.25-0.50 of the stellar mass (JEvans fe Kochaneklll989l ; lAval et al.ll2000l ). 



More massive black holes have Schwarzschild radii exceeding -Rtidaii stars cross the horizon intact 
and / acc = 1. 

For encounters with a binary star having semimajor axis at,i n an d compone nt masses m\ 
and iri2, the probability of a capture/ejection is P< 



cap 



1 - D/175 (JHiUsI Il98fih . where D 



(-R c i ose /abm)/("i,M.), -Rciose is the distance of closest approach, and f(m, M.) 
of M. and mi + m<i . Conservatively setting P cap 
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7ri?^ lose , interact ions with bina r ies (ahin <", 0-1 AU) a re > 100 times more 



Adopting a, 

likely than interactions with single stars (JHopmanl 12003 ; lAntonini et al.ll2011 



The lHilla (J1988I ) proposal yields (i) ejected stars traveling o ut into the Galaxy and ( ii) captured 
stars bound to the black hole. For M. ~ 4x 10 6 M in Sgr A* (JEisenhauer et al.ll2005l ; iGhez et al. 
20081 ). expected ejection speeds are ~ 1000-15 00 km s~ l . Decelera t ion through the Gal axy reduces 
speeds to ~ 500-1000 km s^ 1 at 40-100 kpc (JBromlev et alJbood : iKenvon et alJbood ). Predicted 



orbits for the resulting bound stars have semimajor axes a 
0.95-1.0. 



1000-10000 AU and eccentricities e 



Observations reveal both HVSs and captured stars. Brown et al.l (J2005I ) discovered SDSS 
J090745. 0+024507, a star with a spectral type (B9), age (< 350 Myr), metallicity ([Fe/Hl~ 0) 
radia l velocity (850 km s _1 ), and distance (~ 70 kpc) consistent with Hills' prediction (IBrown et al 
20091 ) . Surveys cove ring roughly 25% of the sky reveal ~ 20 known HVSs with speeds sufficient to 



escape the Galaxy (JEdelmann et ajJ|2005J ; iHirsch et al J 120051 : IBrown et al. 



properties of these stars are consistent with a Galactic Center origin (JBrown et al.l I2O120 



2009 



2012) Observed 



At t he Galactic Center, the S star cluster is plausibly composed of the boun d former partners 
of HVSs JGould fc QuillenlboO.d : iGinsburg fc Loeblbood : b'Learv fc L oeblbood). The ~ 20 ma s- 



sive s tars in this cluster have orbits with a < 80 00 AU and a broad range of e (JEckart &: Genzel 



1997 



Ghez et al. 



1998 



2008 



Genzel et al.ll2010l ). The semimajor axes are consistent with Hills' 



predictions. Although observed ecc entricities are sma ller than predicted, orbital diffusion after 
capture can reduce e substantially (jPerets et al.l |2009| ) . The origin of these stars — or their bi- 
nary progenitors— is uncertain. If they for med in young stellar disks in the Galactic Center 



(JLockmann et al.l 120081 ; iMadigan et al.l |2009| ). their number count and orbital distribution may 
represent an atypical snapshot of the Galactic Center. A clear connection between the S stars and 
HVSs is key to testing the binary disruption scenario. 

Measured produc tion rates for HVSs and S stars suggest a similar formation mechanism (e.g., 
O'Learv & Loebl bood ) ■ Known HVSs have masses of 2.5-3.0 M and travel times of 60-200 Myr 
( Brown et al.ll2012l ). The production rate is roughly 1 x 10 -6 yr _1 . Correcting this rate for unob- 
served, lower mass stars using a standard (differential) initial mass function (IMF) £(m) oc m~ q ~ 1 
with q = 1.00-1.35, the production rate for HVSs of all masses is 2-8 x 10~ 5 yr -1 . S stars have 
masses exceeding 5 M and main sequence lifetimes t ms < 100 Myr (IGhez et al.lll998l : I Genzel et al. 



201Q ). yielding a production rate of 2 x 10 7 yr . Accounting for lower mass stars implies a rate 



of 1-4 x 10~ 5 yr -1 . Within the errors, the rates implied by HVSs and S stars agree. The good 
agreement between the binary disruption rates derived from two populations with very different 
lifetimes suggests a common, steady production mechanism. 

For a more realistic est imate of the bina ry disruption rate from HVSs and S stars, we assume 
continuous star formation (JFiger et al.l 12004 ). and steady, random scattering of binary stars to- 
ward the Galactic Center. Thus the pool of progenitors, accumulated over 10 Gyr, contains more 
older, lower mass stars than predicted by the IMF. To compensate for this population, we adopt 
t ms < 500 Myr for HVSs and t ms < 100 Myr for S stars. The inferred production rates increase by 
~ 10 Gyr/500 Myr = 20 for HVSs and ~ 10 Gyr /100 Myr = 100 for S stars. Thus, the maximum 
rate for binary disruptions in the Galactic Center is roughly 1-2 x 10~ 3 yr -1 . 

This range of rates agrees with predictions. In simple models, binaries interact with the central 
black hole at a rate of n (i>.f g <7,), where n is the number density of binaries, v ~ 100 km s _1 is 



a characteristic speed ( Figer et al.ll2003l ). and / g cr, ~ lx 10 8 AU is the cross-sectional area for 
-^ciose ~ 100 AU (eq. [2]) and a gravitational focusing factor f g . For a central mass density po, 
n ~ fbinPo / fn&vg (1 + /bin); where /bin is the fraction of stars in short period binaries and m avg sa 
0.3 M is the average mass of a star. In the solar neighborhood, the fraction of stars in binaries 
of all periods is /bin = 0.5-0.7 (lAbtlll983r). For the short period binaries likely to pro duce HVSs 
(P < 1 yr), /bin ~ 0.1 (JYu &; Tremaindl2003l ). The total disruption rate is then (see also iHills 
Yu fc Tremainelhoosh 



1988 



^dis 



/bin 



P0 



1 "r* /bin ^-avg 
10 -3 yr -l 



(vfgcr.) 



for po = 1.4 x 10 4 M© pc~ 3 (e.g., Bromley et alJl200d ) and / bin 
maximum rate of binary disruption implied by observations. 



(3) 

(4) 
0.1. This rate agrees with the 



The maximum disruption rate assumes that scattering processes rapidly fill the "loss cone," 
the subset of Galactic Center orbits that pass close to the black hole. If scattering supplies binaries 
to the loss cone more slowly than disruption removes them, the di sruption rate is compar able to 
the scattering rate. When two-body relaxation fills the loss cone, lYu fc Tremaind (|2003l ) derive 
fcdis ~ 10~ 5 yr -1 for /hin = 0-1. Because oth e r processes can fill the loss cone, this estimat e 



is probably a lower limit (JMerritt &: Poon 



2004 



Hopman fc Alexander! 120061; 



When the gravitational potential is nonaxisymmetric and orbits are chaotic, iMerritt fe Poonl ((2004) 



Perets et al. 



2007). 



derive factor of ten larger scattering rates. The expected disruption rate is then k, 
Taken together, the broad range of theoretical estimates, fedis 
the binary disruption rates implied by observations. 



dis 



10" 4 yr" 



10 5 -10 3 yr , is consistent with 



3. Long-term Orbital Evolution of the S Stars and Tidal Disruption Events 



The inferred binary disruption rate agrees with estimates for TDEs from nearby galaxies. 
Observations of two candidate TDEs from long-term ob servations of galaxies in the Sloan Digital 

Candidates from 



Sky Survey suggest one TDE per galaxy ev ery ~ 10 5 yr ( van Velzen et al. 



GALEX data indicate a similar frequency (JGezari et al. 



2011 



20091 ). The good agreement between the 



rates of binary disruption and stellar disruption plausibly suggests a common origin. To explore 
this connection, we consider the evolution of the orbits of bound stars. 

Over 10 Gyr, binary disruption places ~ 10 5 -10 7 bound stars in orbit around the central black 
hole. Just after capture, each star has a semimajor axis and eccentricity 
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for M, = 4 x 10 6 Mq. These equations follow from conservation of energ y and angu l ar momentum, 



with coefficients from numerical simulations of binary disruptions (e.g., lHillslll988l : iBromley et al 
200fil ). 



Gravitational interactions among the bound stars produce changes in Obnd and Cbnd on a 
relaxati on time T rp ] . Stars w i th t m .i S, t vp \ explode as supernovae or evolve into red giants and white 



dwarfs (Frank fc Rees 1976; Hils fc Bender 



1995 



Hopman fe Alexander 



2006 



Merrittil20ld ). Less 



massive stars may evolve onto orbits w ith periastr on distances a per i < -Rtidai- These stars produce 
TDEs and accrete onto the black hole (|Reeslll988l ). 



The rate of TDEs depends on r re i- If two-body scattering dominates, the relaxation time is 
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when ~ 10 6 stars with an average mass of 0.3 M orbit at ~ 5000 AU with orbital period 
P. Other processes, such as resonant relaxation, probably shorten the relaxation time (e.g., 
Hopman fc Alexander! 120061 ). Captured stars wander through < e < 1 on this time scale. 



The relaxation time sets two properties for the bound population. Stars with t ms < r re i 
(m > 2-4Mq) evolve into supernovae or white dwarfs before they encounter the black hole. For a 
normal IMF, the fraction of stars that suffer this fate is a few per cent. The rest of the population 
establishes an approximate steady-state, where the black hole accretes a star for every captured 
star. For the current mass of Sgr A*, this process requires ~ 0.1-1 Gyr. After this period, there is 
a tidal disruption event every 10 3 -10 5 yr. 



Numerical simulations do not consider the long-term evolution of ~ 10 5 -10 7 low mass stars 
bound to a ce ntral black hole. H owever, smaller simulations confirm the general features of this 



evolution (e.g.. iPerets et al.ll2009l ). Gravitational interactions lead to global changes in a and e on 
the relaxation time r re i. At least 20% to 30% of bound stars suffer TDEs. Thus, Sgr A* should 
accrete a 0.1-0.5 M star every 10 3 -10 5 yr. 



4. Growth of the Central Black Hole 

The upper limit on the tidal disruption rate can add mass comparable to the total mass of 
Sgr A* . To explore how binary disruption impacts M, , we consider a model where the capture rate 
is the collision rate in equation ([3]). For a constant central density po = 10 4 M pc~ 3 , velocity 
v = 100 km s _1 , and binary fraction /bin = 0.1, our time-dependent calculations begin with an 
initial M. = M see( j. Captured stars accumulate until time t = r Te \ and then accrete with efficiency 
/ acc onto the black hole. The growth of M. depends only on the initial mass and the accretion 
efficiency. 

Figure [U shows results for several seed masses and / acc = 0.25 and 1.00. In FigureHJ small seed 
masses and / acc < 1 yield modest growth rates over 10 Gyr. Larger seed masses of a few x 10 5 M 
yield M. w 4 x 10 6 M within 10 Gyr. 

To show how black hole growth depends on po and /bin, we define a growth parameter 

" = (¥)(£)(&) (9) 

where p c is the central mass density near the black hole. Thus, rj = 1 corresponds to models with 
the Milky Way's central stellar mass density. Calculations with rj > 1 (rj < 1) yield larger (smaller) 
growth rates. 

Figure [2] displays M, at 10 Gyr as a function of r\ and the initial seed mass. In this model, 
massive black holes require large accretion efficiency, large central density, and a large fraction of 
short period binaries. Low values for r/ yield low M,. Figs. [TH2] imply that the current binary 
disruption rate implies a substantially lower mass for Sgr A* in the past. 

To explore black hole growth in more detail, we add the evolution of i?tidai (eq. [I] ) and -R c iose 
(eq. [2]) as the black hole grows in mass. Large gravitational focusing factors imply that the 
disruption rate is independent of v (eq. [3]). We also allow the central density to vary i n time, 



adopting p(t) oc M, 7 and 7 = 1-2, as observed in elliptical galaxies (JFaber et al.lll997l ). This 
density is bounded by a maximum initial density of p max = 10 6 M pc -3 . 

In this model, the mass accretion rate is 

dM, _ 2/ acc m avg (1 + e tida i) 

dl T re l 



which includes the rate of diffusion into orbits that exceed e t idai) the eccentricity corresponding to 
a perihelion distance of -Rtidai- The resulting growth of the black hole is not sensitive to the exact 
form of this equation. As in the simple model (Figs. dH2|)) the rate invariably adjusts to achieve a 
steady state with a balance between capture of stars and the tidal shredding by the black hole. 

Although there may be no direct causal connection between the black hole and the central 
density outside of its sphere of influence, our approach implies that the central density diminishes 
in time as M. grows. Once the central density falls below p max , the numerical results suggest that 
the central density declines with time roughly as p(t) ~ t~@, where /3 ~ 1. 

Figure [3] summarizes our results. For an initial black hole seed mass of 10 5 Mq, the two curves 
show the central density as a function of black hole mass at t = 10 Gyr. For any M, at 10 Gyr, 
models with 7 = 1 require larger central stellar mass densities than models with 7 = 2. Observation s 



show considerable scatter about the general trend derived from local galaxies (jFaber et al.lll997l ). 
Although we cannot discriminate among plausible sets of variables, the growth models yield a fair 
representation of the trend among galaxies with large central stellar mass densities. These models 



also appear to fare better than the power-l aw correlation p oc M. from dynamical arguments 



(dashed line in the figure; iFaber et al.lll997l ). 



5. Summary 

We develop a model for binary disruption in the Galactic Center. The model relates observa- 
tions of HVSs in the halo and the S star cluster in the Galactic Center, TDEs in nearby galaxies, 
and the growth of central black holes. The black hole in the Galactic Center disrupts close binaries 
at distances < i? c iose (eq- 0). One companion (the HVS) is ejected at high speeds; the other 
(S star) becomes bound to a black hole. Among the bound stars, dynamical interactions push 
some stars inside the tidal radius (eq. [1]). Shredding of stars inside -Rtidai produces TDEs. After 
a relaxation time (eq. [7]), the rate of tidal disruptions roughly equals the capture rate. For large 
enough capture rates, binary disruptions add significant mass to the central black hole. 

Current data suggest this picture is plausible. Observations of HVSs and S stars yield consistent 
production rates of 10~ 5 -10~ 3 yr _1 ; these rates agree with theoretical estimates. For t vc \ ~ 0.1- 
1 Gyr, the expected rate for TDEs, 10 -5 -10 -3 yr , is close to current estimates from small samples. 
More plentiful singl e stars may contribute to TDEs but the rate is at or below the level of the binary 



disruption channel (jBrockamp et al.ll2011l ). For rates associated with binary disruption, Sgr A* has 
grown by at least a factor of 2-4 in the past 5-10 Gyr. At higher rates, binary disruption can 
produce black holes with M. > 10 8 Mq. This growth model yields clear relations between the final 
M. and the central stellar mass density. These relations are reasonably consistent with existing 
data. 

This model is attractive in connecting HVS and S stars in the Milky Way to TDEs in nearby 
galaxies. To test this connection, A-body simulations with larger numbers of stars bound to the 



black hole would improve estimates of the relaxation time and relations between the capture rate 
and the disruption rate. Theoretical investigations of TDEs involving the bound former partners 
of disrupted binaries might yield different light curves which could be tested by observations. 

Several observations also provide essential tests. 

• Extending HVS surveys to the southern sky with SkyMapper would improve estimates of the 
HVS production rate. In the future, LSST should detect candidates with lower mass and at 
larger distances in the halo. Spectroscopic confirmation with current 6-10-m telescopes and 
proposed 25-50-m telescopes would yield even better constraints on the HVS population. 

• Discovery of lower mass stars among the bound star population yields another test of the 
binary disruption hypothesis. Deeper observations of the Galactic Center with existing large 
telescopes, the James Webb Space Telescope, and planned 25-50-m telescopes should reveal 
this population. 

• Larger samples of candidate TDEs enable better comparisons of the binary disruption rate 
in the Milky Way with tidal disruption rates in external galaxies. Aside from existing sur- 
veys with the SDSS and Pan-Starrs, surveys with SkyMapper and LSST will yield strong 
constraints. In any picture, low mass stars should produce most TDEs. With constraints on 
the black hole mass, deriving a luminosity function for this process might place estimates on 
the mass function of disrupted stars. Rates as a function of central stellar mass density tests 
theories for the capture rate. 

Calculations of black hole growth by tidal disruption predict a correlation between the stellar 
mass (velocity dispersion) and the final M.. Although our models cannot predict the observed 
M. oc a relation, more sophisticated calculations could reveal M.-ct correlations which might be 
tested observationally. 
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Fig. 1. — Growth of a black hole from a seed mass. Curves correspond to seed masses between 
10 5 M and 3.5 x 10 6 M . The lower (upper) curves have / acc = 0.25 (/ acc = 1). 
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Fig. 2. — Tracks of constant final M, as a function of the growth parameter r\ and the initial M,. 
For the accretion rates expected in the Galactic Center, initial masses of ~ 10 6 M yield the mass 
of Sgr A* in 10 Gyr. 



13 



10 5 



10 4 



O 



© 1000 



o 



100 



10 



1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



10 




J I I 



10 8 10 9 

M BH (M ) 



Fig. 3.—- Central stellar mass density as a function of final black hole mass. For models with a 
seed mass of 10 5 M , / acc = 0.25, and p{t) oc M.~ 7 , the upper (7 = 1), middle (7 = 1.5), and lower 
(7 = 2) solid lines show predicted relations between po and M. at 10 Gyr. We generate each curve 
by varying the initial density p(t = 0), and letting p(t) evolve as M, grows. Thus, while a large 
initial density yields a large black hole mass, the present-day density po is drive n to a small value i f 
7 > 0. Observations s upport this correlat i on. F illed circles show data from the lFaber et al.l (J1997I ) 
sample with M, from lAller fc Richstond (|2007l ) . The open circles are from the sa me sample, but 
with the assumption that M, ~ -i> 4 , where v is the measured velocity dispersion ( Tremaine et al. 
2002). The dashed line is a power-law relation with 7 = 1.3, representative of galaxies in the local 
universe. 



